AULA 


Modelos matemáticos 
de E.D.O. de primeira ordem 


META: 

Descrever, resolver e interpretar alguns modelos matemáticos de 
fenômenos reais que utilizam as Equações Diferenciais Ordinárias 
de primeira ordem. 

OBJETIVOS: 

Ao fim da aula os alunos deverão ser capazes de: 

Modelar matematicamente alguns fenômenos que utilizam Equações 
Diferenciais Ordinárias de primeira ordem. 

Interpretar as soluções dessas equações. 

PRÉ-REQUISITOS 


Os conhecimentos das quatro aulas anteriores. 


Modelos matemáticos de E.D.O. de primeira ordem 
5.1 Introdução 


Caros alunos nessa quinta aula aplicaremos os métodos de re- 
solução de E.D.O.'s aprendidos nas aulas anteriores e assim vere- 


mos como essas equações de fato estão presentes na nossa realidade. 


5.2 Dinâmica populacional 


Em se tratando de dinâmica populacional existem dois modelos 
matemáticos, muito bem sucedidos, usados para descrever o cresci- 
mento ou decaimento de uma população. 

O primeiro modelo, proposto em 1798, é devido ao economista in- 
glês Thomas Malthus, conhecido como modelo Malthusiano. Nesse 
modelo, Malthus se baseia na hipótese de que a taxa de variação de 
uma população em um determinado instante é proporcional à po- 
pulação total naquele instante, ou seja, se tomarmos as variáveis t 
como tempo e P(t) como a quantidade de indivíduos da população 


no instante t, matematicamente, Malthus quer dizer que 


onde dP/dt é a taxa de variação da população e k é a constante 
de proporcionalidade. Esse modelo se mostra eficaz apenas para 
estudar pequenas populações em curto espaço de tempo, como por 
exemplo, crescimento de micro-organismos que se reproduzem por 
mitose. 

O segundo modelo, proposto em 1838, é devido Verhulst, matemático 
belgo. Esse modelo não considera k como uma constante uma 
vez que o crescimento de uma população aciona automaticamente 


certos mecanismos de controle visando reduzir a taxa de cresci- 


( so b 
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mento. O modelo proposto por Verhulst consiste em supor que 
k decresce linearmente com o aumento da população, ou seja, 
k=a-bPa,b e R. Assim, segundo Verhulst, a dinâmica de 
uma população pode ser estudada pela equação 


dP 
-— = (a-bP)P. 
E (a — bP) 


Esse ainda não é o modelo que considere todas as situações reais 
possíveis, contudo seus resultados são melhores do que os do mo- 
delo de Malthus. 

Observe que tanto o modelo de Malthus quanto o modelo de Ver- 
hulst são do tipo que se resolvem pelo método de equações sepa- 
ráveis. Vamos aos cálculos !! 


Analisemos primeiramente o modelo de Malthus 
— = kP(t) sujeito à condição P(to) = Po. (5.30) 


Separando as variáveis em (5.30), obtemos 


dP 
= kdt, P(to) = Po. 


E integrando, segue que P(t) = Pye*tt-to). Observe que nesse mod- 
elo o crescimento (k > 0) ou decrescimento (k<0) de uma popu- 
lação é exponencial! 

Consideremos agora o modelo de Verhulst 


dP 
E (a— bP)P, sujeito à condição P(to) = Po. (5.31) 


Usando frações parciais, temos que 


1 1 b 


(a-bP)JP aP' a(a-bP) 
assim o P.VI. (5.31) se reduz a 
P P 
É d = di Po Dj 


aP ii a(a — bP) 
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Integrando a E.D.O. acima, temos 


1 1 
In|P| inlja-bP|=t+c,R. 
a a 


Daí, se P(to) = Po, fazendo algumas simplificações, chegamos a 


2" aPg 
“bP+4(a-bPoe-t-to)' 


P(t) 


Abaixo descrevemos a curva integral da solução 


Exemplo 5.1. A população de uma cidade cresce a uma taxa 
proporcional à população presente num instante de tempo t. À 
população inicial de 500 cresce 15% em 10 anos. Qual será a pop- 
ulação em 30 anos? 

Pelo enunciado vemos que esse problema de crescimento popula- 
cional leva em conta o modelo de Malthus, assim a expressão da 
solução é dada por 


P(t) = cef',ce R. 


Sabemos, pelo enunciado, que P(0) = 500 e P(10) = 500 + 0,15- 
500 = 500 - 1,15. Assim, 500 = ce? que resulta c = 500 e 500 - 
1,15 — 500e!º* que nos dá k = nu. Assim, a população em 30 


e 760. 


ln 1,15 
0 


anos será P(30) = 500e% "1 


5.3  Datação da idade de um fóssil 


A datação de um fóssil usando carbono, desenvolvida pelo químico 
Willard Libby, por volta de 1950, baseia-se no fato de que o isó- 
topo carbono 14 é produzido na atmosfera pela ação da radiação 
cósmica sobre o nitrogênio. A razão da quantidade de C— 14 em re- 
lação ao carbono comum na atmosfera parece ser uma constante e, 


consequentemente, a quantidade proporcional de isótopo presente 


(82, 
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em todos os organismos vivos é a mesma na atmosfera. Quando 
um organismo morre ele deixa de absorver C — 14, dessa maneira 
comparando a quantidade proporcional de C — 14 presente em um 
fóssil com a razão constante encontrada na atmosfera e sabendo 
que o tempo de meia vida (tempo necessário para que a metade 
dos átomos de uma quantidade inicial de uma substância se desin- 
tegrem ou transformem-se em átomos de um outro elemento) do 
radiativo C — 14 é aproximadamente 5600 anos, é possível estimar 
o ano que aquele organismo vivo morreu, (ZILL, 2003). 

Observe, pelo exposto acima, que o problema da datação de um 
fóssil é basicamente estruturado pelo dacaimento radioativo do isó- 
topo C— 14 presente nos organismos vivos. Acredita-se que o fenô- 
meno do decaimento radioativo acontece de maneira que a taxa de 
decaimento segundo a qual o núcleo de uma substância decai é 
proporcional à quantidade da substância remanescente no instante 


t, ou seja, matematicamente falando 


onde A(t) é a quantidade de substância remanescente no instante 
de tempo t. Dessa maneira obtemos uma equação diferencial semel- 
hante ao do modelo populacional de Malthus. Na próxima seção 
veremos um outro fenômeno que é modelado matematicamente 


pelo mesmo tipo de E.D.O.. 


Exemplo 5.2. No sertão paraibano foram encontradas várias pin- 


turas gravadas em pedra. Arqueologistas usaram pedaços de madeira 


queimados nesse local para datar a idade dessas pinturas pré- 
históricas. Com a ajuda de um aparelho de medição, detectou-se 


que 90% do € — 14 havia decaído daqueles pedaços de madeira 
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encontrados. Estime a idade das pinturas encontradas. 


5.4 Esfriamento e aquecimento de um corpo 


De acordo com a lei do aquecimento e resfriamento de um corpo, 
descrita por Newton, a taxa segundo a qual a temperatura de um 
corpo varia é porporcional à diferença entre a temperatura do corpo 
e a temperatura do meio que o rodeia (temperatura ambiente). 
Assim, considerando T(t) a temperatura do corpo no instante de 
tempo t e Ty a temperatura do ambiente em que se encontra o 
corpo, temos a seguinte equação que modela matematicamente o 


fenômeno do aquecimento ou resfriamento de um corpo 
— =k(T(t)-T),kER. 


Convidamos o aluno interessdo a pensar no seguinte problema. 


Mãos a obra: Um indivíduo é encontrado morto em seu escritório 
pela sua secretária que liga imediatamente para a polícia. Quando 
a polícia chegou imediatamente mediu a temperatura do defunto, 
achando 35ºC; uma hora depois mediu novamente obtendo 34, 2ºC. 
Sabendo que a temperatura do escritório é de 20ºC€ e supondo que 
a temperatura normal de uma pessoa viva seja constante e igual a 


36,5º0, 


a- Determine quanto tempo depois da morte do indivíduo a polícia 


chegou no escritório. 


b- Sabendo que a polícia chegou 2 horas depois que a secretária 


telefonou, decida se o indivíduo morreu antes da secretária tele- 
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fonar para a polícia ou depois da secretária ter telefonado para a 


polícia? 


5.5 Circuitos elétricos 


Considere um circuito em série de malha simples, descrito abaixo, 
contendo um resistor com resistência R, um indutor com indutân- 


cia L e um capacitor com capacitância C. 


Figura 5.1: Circuito em malha simples. 


Denotemos por i(t) a corrente que passa pelo circuito quando a 
chave é fechada e q(t) a carga no capacitor em um instante t. De 
acordo com a Lei de Kirchhoff, a voltagem aplicada E(t) em 
uma malha fechada deve ser igual à soma das quedas de voltagem 
na malha. Sabe-se que quando uma corrente i(t) passa por um 
resitor ou um indutor ou um capacitor, a tensão retida por esses 
componentes é dada, respectivamente, por R-1, L- E e õ “q. Dessa 


maneira, usando a Lei de Kirchhoff, obtemos 


dá 1 
psi: 
(6) rr É 


s d Es : z , 
e como 1 = = esta equação pode ser reescrita da seguinte maneira 


d dg 1 
SALG toa (5.32) 


E ai Era 
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Portanto, a quantidade de carga ao longo do tempo num circuito 
de malha simples, como o descrito acima, é dada pela solução da 


E.D.O. (5.32). 


Mãos a obra Uma força eletromotriz 


0, t > 20 


é aplicada em um circuito em série RL no qual a indutância é 20 


henrys e a resistência é de 2 ohms. Ache a corrente i(t) se i(0) = 0. 


5.6  Diluição de soluções 


Suponha que, em um tanque se tenha 50 litros de água pura, a 
qual começa a receber uma salmoura (solução de água com sal) de 
concentração 3 kg de sal por litro de solução a uma razão constante 
de 2 litros por segundo. Um mecanismo de agitação no tanque 
mantém homogênea a solução que vai sendo formada. Começa a 
se retirar do tanque a solução formada, na razão de 2 litros por 
segundo. Como podemos saber a quantidade de sal presente no 


tanque no instante de tempo t? 


Observação 5.1. A mistura de duas soluções salinas com concen- 
trações diferentes também é feita usando a mesma idéia desse tipo 


de questão. 


Seja x(t) a quantidade de sal em kg presente no tanque num tempo 
t. Portanto a concentração de sal na solução é 59kg/1. Daí tem-se 


que a taxa de variação de sal no tanque obedece a seguinte E.D.O. 


dx x 
e S=3 Ro (5.33) 
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onde TE e TS quer dizer, respectivamente, taxa de entrada de sal 
e taxa de saída de sal. 
É fácil ver que a equação (5.33) é do tipo linear ou se preferir 


podemos resolvê-la separando as variáveis. 


5.7 Trajetórias ortogonais 


Antes de falarmos como obter trajetórias ortogonais é necessário 


que definamos alguns termos que usaremos aqui. 


Definição 5.1. (Curvas Ortogonais)- Duas curvas C, e C» são 
ortogonais em um ponto P se , e somente se, suas retas tangentes 


T1 e T3) forem perpendiculares no ponto de interseção, ou seja, se 
Memo — —1, (5.34) 


onde my e ma são respectivamente os coeficientes angulares das 
retas tangentes. (aqui estamos supondo retas tangentes não pa- 


ralelas aos eixos coordenados.) 


Definição 5.2. (Trajetórias ortogonais)- Quando todas as cur- 
vas de uma família G(x,y,c1) = O interceptam ortogonalmente to- 
das as curvas de uma outra família H(x,y,c>) = 0, dizemos que 


as famílias são trajetórias ortogonais uma da outra. 


Exemplo 5.3. A família de retas y = c;x e a família de círculos 
concêntricos «2 + y? = c» são trajetórias ortogonais. Veja figura 


abaixo 


Dada uma família de curvas G(x, y,c1) = O para encontrarmos uma 


família de curvas H(x,y,co) = O tal que G e H sejam trajetórias 
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Figura 5.2: Trajetórias ortogonais. 


ortogonais devemos proceder da seguinte maneira: suponha que 


G(x,y,c1) = O seja a família de soluções da E.D.O. 


a = g(x,4), 

então as inclinações das retas tangentes às soluções dessa família no 
ponto (x,y) são dados por g(x,y). Portanto, segundo a condição 
(5.34), para acharmos a família de curvas ortogonais as curvas 


G(x,y,c1) = 0, devemos achar a família de soluções da E.D.O. 


1 
g(w,y)* 


HE 
3 


Exemplo 5.4. A família de hipérboles y — “ é uma família de 


soluções da E.D.O. 
dy y 


de q 
Assim, a E.D.O. da família ortogonal a essa é dada por 


dy x 
dr y 
A família de soluções dessa E.D.O. é dada por y? — «2 = c» que 
também é uma família de hipérboles com eixo real sobre o eixo y. 


Abaixo descrevemos as duas famílias 


Famílias de curvas mutuamente ortogonais aparecem de modo na- 


tural em muitos problemas físicos. Por exemplo, as curvas de fluxo 
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Figura 5.3: Trajetórias ortogonais. 


de calor numa lâmina fina de metal são ortogonais às curvas de 
temperatura constante, as chamadas curvas isotérmicas. 

Ainda existem mais aplicações a serem estudadas e exploradas. 
Para quem deseja conhecer mais aplicações e se aprofundar nesse 
assunto, sugerimos o livro Equações Diferenciais com aplicações de 


Bassanezi e Ferreira Jr., Editora Habra. 


5.8 Conclusão 


Na aula de hoje, vimos diversas aplicações para as E.D.QO.'s de 
primeira ordem. É muito amplo o campo de aplicabilidade dessas 
equações, as quais são usadas em problemas de Biologia, Física, 
Química e Matemática. Alguns problemas vistos nessa aula são 
modelados, basicamente, pela mesmo tipo de equação. É o caso 
do modelo de crescimento populacional Malthusiano, esfriamento 


ou aquecimento de um corpo e datação da idade de um fóssil. 
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RESUMO 


Hoje vimos diversos problemas práticos, tais como: 


1. Estudo da dinâmica populacional 

2. Datação da idade deum fóssil 

3. Esfriamento ou aquecimento de um corpo 
4. Circuitos elétricos 

5. Diluição de soluções 


6. Como obter trajetórias ortogonais. 


Todos esses problemas são modelados por equações diferenciais 
ordinárias de primeira ordem cuja resolução aprendemos nas aulas 


03 e 04 desse livro. 


PRÓXIMA AULA 


Em nossa próxima aula começaremos expondo alguns resultados 
teóricos que nos serão imprescendíveis na hora de resolvermos al- 


gmas E.D.O. lineares de ordem maior do que um. 


ATIVIDADES 


Atividade. 5.1. Uma substância radioativa decompõe-se a uma 
razão proporcional à quantidade presente, e no fim de 1500 anos 
reduz-se à metade da quantidade original. Em quantos anos, a 
quantidade original se reduz a três quartos? Qual a quantidade da 


substância encontrada no fim de 2000 anos? 
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Atividade. 5.2. Encontre as trajetórias ortogonais a z+y = cy e” 


que passam por (0,5). 


Atividade. 5.3. Sabe-se que a população de uma certa comu- 
nidade cresce a uma taxa proporcional ao número de pessoas em 
qualquer instante. Se a população duplicou em 5 anos, quando ela 


triplicará? Quando ela quadruplicará? 


Atividade. 5.4. O isótopo radioativo de chumbo, Pb-209, decai 
a uma taxa proporcional à quantidade presente no instante t e tem 
uma meia-vida de 3,3 horas. Se houver 1 grama de chumbo ini- 


cialmente, quanto tempo levará para que 90% do chumbo decaia? 


Atividade. 5.5. Uma força eletromotriz é aplicada a um circuito 
em série LR no qual a indutância é de 0,1 henry e a resistência é de 
50 ohms. Ache a curva i(t) se a corrente inicial é zero. Determine 


a corrente quando t — oo. 


Atividade. 5.6. Deduza que a equação diferencial governando a 
velocidade v de uma massa em queda sujeita à resistência do ar 


proporcional à velocidade instantânea, é 


onde k > 0 é uma constante de proporcionalidade, considere o 


sentido positivo para baixo. Use essa equação para responder 
a)Resolva a equação sujeita à condição inicial v(0) = vo. 


b) Use a solução do item a) para determinar a velocidade-limite 


ou terminal da massa. 


Atividade. 5.7. Um termômetro é levado para fora de um quarto, 


onde a temperatura ambiente é de 5ºF . Após 1 minuto, o ter- 
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mômetro marca 55ºF, e após 5 minutos, 30ºF. Qual era a temper- 


atura inicial no quarto? 
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